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A triangulation of a planar point set S is a maximal plane straight-line graph with vertex set S. In 
the minimum-weight triangulation (MWT) problem, we are looking for a triangulation of a given 
point set that minimizes the sum of the edge lengths. We prove that the decision version of this 
problem is NP-hard, using a reduction from PLANAR 1-IN-3-SAT. The correct working of the 
gadgets is established with computer assistance, using dynamic programming on polygonal faces, 
as well as the /3-skeleton heuristic to certify that certain edges belong to the minimum-weight 
triangulation. 
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1. INTRODUCTION 

Given a set S of points in the Euclidean plane, a triangulation T of S is a maximal 
plane straight-line graph with vertex set S. The weight of T is defined as the total 
Euclidean length of all edges in T. A triangulation that achieves the minimum 
weight is called a minimum-weight triangulation (MWT) of S. 

The problem of computing a triangulation for a given planar point set arises nat- 
urally in many applications such as stock cutting, finite element analysis, terrain 
modeling, and numerical approximation. The minimum-weight triangulation has 
attracted the attention of many researchers, mainly due to its natural definition of 
optimality, and not so much because it would be important for the mentioned ap- 
plications. We show that computing a minimum-weight triangulation is NP-hard. 
Note that it is not known whether the MWT problem is in NP, because it is an open 
problem whether sums of radicals (namely, Euclidean distances) can be compared 
in polynomial time [Blomer 1991 J . This problem is common to many geometric 
optimization problems, like, perhaps most famously, the Euclidean Traveling Sales- 
person Problem. To get a variant of the problem that is in NP, one can take the 
weight of an edge e as the rounded value HMhl- Our proof also shows that this 
variant is NP-complete. 

Our proof uses a polynomial time reduction from POSITIVE PLANAR l-IN-3- 
SAT, a variant of the well-known PLANAR 3-SAT problem, which is a standard 
tool for showing NP-hardness of geometric problems. 



1.1 Optimal Triangulations 

Usually, a planar point set has (exponentially) many different triangulations, and 
many applications call for triangulations with certain good properties (see Figure H]). 
Optimal triangulations under various optimality criteria were extensively surveyed 
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Fig. 1. A planar point set and different ways to triangulate it. The greedy triangulation (a) 
is constructed incrementally, always adding the shortest possible edge. In this example, it is 
shorter than the Delaunay triangulation (b), which avoids skinny triangles. A minimum weight 
triangulation for the point set is shown in (c). The length L of the MWT can decrease when 
additional Steiner points are allowed (d). 



by Bern and Eppstein 1992 



The Delaunay triangulation is perhaps the best known triangulation. It arises as 
the dual of the Voronoi diagram and can be computed optimally in O(nlogn) time, 
using classical techniques |de Berg et al. 2000| . Since it simultaneously optimizes 
many objective functions (it maximizes the smallest angle, it minimizes the maxi- 
mum circumcircle, as well as the maximum smallest enclosing circle of all triangles), 
it is often the triangulation of choice. 

Certain other criteria, such as the minimum maximum angle, the maximum min- 
imum height, or the minimum maximum distance of a triangle from its circum- 
center, can be optimized in polynomial time using the edge insertion technique 
[Bern et al. 1 993; Ed elsbrunner et al. 1992j . a sophisticated application of dynamic 
programming. 

It makes a difference whether it is allowed to add new points, so-called Steiner 
points, to the planar point set. As Figure Q] shows, this can sometimes help to 



improve the objective function. Bern and Eppstein [1992 give many similar ex- 



amples. Here we only mention one result by Eppstein [1994| , who showed that the 
minimum weight Steiner triangulation can be approximated up to a constant factor 
in 0(n log n) time, using quadtrees. 

1.2 History of the MWT problem 

The minimum-weight triangulation problem dates back to the 1970s and has been 
called "perhaps the most longstanding open problem in computational geometry" 



|Bern and Eppstein 1992] , At the end of the classical book of Garey and Johnson 1979] 
on NP-completeness, there is a list of 12 major problems whose complexity status 
was open at the time of writing. With the present paper, ten problems from this list 
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have been resolved by proving NP-hardness or by exhibiting a polynomial-time al- 
gorithm (see [Johnso n 2005] for a recent status update on the list). As of now, only 
two problems from the original list remain open, namely Precedence-Constrained 
3-Processor Scheduling and the notorious Graph Isomorphism problem. 



Early attempts. It seems that the MWT problem was first considered by Diippe and Gottschalk [1970 



who proposed a greedy algorithm which always adds the shortest possible edge to 
the triangulation. Later, Shamos and Hoey [1975] suggested using the Delaunay 



triangulation as a minimum- weight triangulation. |Lloyd 1977] provided examples 
which show that both proposed algorithms usually do not compute the MWT (Fig- 
ure [[]) . He also shows that it is NP-complete to decide whether the edge set of 
a given planar straight-line graph (with crossing edges) contains a triangulation. 
After this, countless researchers attacked the MWT problem from many different 
angles. In one line of attack, researchers used classical optimization techniques such 
as dynamic programming or branch and bound, but this soon became infeasible. 
In other lines of research, people looked at relaxed variants of the problem: Maybe 
there exist reasonable restrictions of the problem for which efficient algorithms can 
be found, and if the problem cannot be solved exactly, maybe good approxima- 
tions can be computed efficiently. Finally, there was an attempt to gain a better 
understanding of the geometric properties of the MWT in order to find footholds 
for effective heuristics. We now describe these approaches in more detail. 



Dynamic Programming. Gilbert 1979] and Klincsek 1980 1 independently showed 
how to compute a minimum- weight triangulation of a simple polygon in 0(n 3 ) 
time by dynamic programming. In fact, this problem has become one of the stan- 
dard textbook examples (or exercises) for illustrating the dynamic programming 
paradigm. There have also been attempts to attack the general problem with dy- 



namic programming techniques. For example, Cheng et al. [1995 used dynamic 



programming in order to compute a minimum-weight triangulation of a given point 
set S in 0(n k+2 ) time if a subgraph of a MWT of S with k connected components is 



known. Using branch and cut, Kyoda et al. [1997 managed to compute MWTs of 



100 points, but for large point sets mere dynamic programming becomes absolutely 
infeasible. 

Restricted Instances. For restricted classes of point sets, it is possible to compute 



the MWT in polynomial time. For example, Anagnostou and Corneil [1993 gave 



an algorithm to compute the MWT of the vertex set of k nested convex polygons 
in 0(n 3fc+1 ) time. More recently, Hoffmann and Okamoto 2006] showed how to 



obtain the MWT of a point set with k inner points in 0(6 k n b logn) time 

Approximations. In another line of attack, researchers were looking for triangula- 
tions that approximate the MWT. The Delaunay triangulation is not a good candi- 
date, since it may be longer by a factor of Q(n) |Kirkpatrick 1980; Ma nacher and Zobrist 19 79] . 
The greedy triangulation approximates the MWT by a factor of Q{^Jn) [Manacher and Zobrist 1979; 



Levcopoulos 1987 ; Levcopoulos and Krznaric 1996| . Plaisted and Hong 1987] showed 
how to approximate the MWT up to a factor of O(logn) in 0(n 2 log n) time. 



Levcopoulos and Krznaric [1996 1 introduced quasi-greedy triangulations, which ap- 



proximate the MWT within a constant factor. Remy and Steger [2006 discovered 
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an approximation scheme for MWT that runs in quasi-polynomial time: for every 
fixed e, it finds a (1 + e)-approximation in n ( log n ^ time. 

Subgraphs and supergraphs. A different line of research tried to identify criteria to 
include or exclude certain edges of the MWT. iGilbert [1979] showed that the MWT 



always contains the shortest edge. Yang et al. [1994] extended this result by proving 



that edges which join mutual nearest neighbors are in the MWT. A larger subgraph 
of the MWT, the (3-skeleton, was discovered by Keil 19S 
tion 14. 1[ In practice, the /^-skeleton has many connected components 



We describe it in Sec- 
and thus 



does not help much in computing a MWT. Often, the LMT-skeleton heuristic de- 



scribed by Dickerson et al. [1997] yields much better results. It uses the simple fact 



that a MWT is locally optimal in the sense that it cannot be improved by flipping 
the diagonals of a convex empty quadrilateral in the point set. The LMT-skeleton 
made it feasible to compute the MWT for larger, well-behaved point sets, and it has 
been the subject of numerous further investigations [Beirouti and Snoeyink 1998 



Cheng et al. 1996| lAichholzer et al. 1999UBelleville et al. 19961 IBose et al. 2002j . 

Approaching the problem from the other direction, Das and Joseph [1989 de- 
fined the diamond test, which yields a supergraph of the MWT: An edge e can only 
be contained in the MWT if at least one of the two isosceles triangles with base e and 
base angles ir/8 is empty. This constant was improved to 7r/4.6 Drysda le et al. 2001| 
(see Figure [9] below). The diamond test gives an easy criterion to exclude impos- 
sible edges from the MWT. Usually, this eliminates all edges except a set of 0(n) 
remaining candidate edges. (This statement is true for random point sets, with 
high probability. With bucketing techniques, such a set of 0(n) edges can be found 
in linear expected time |Drysdale et al. 19 95 1.) 

1.3 Our Methods and Results 

1.3.1 Historical Perspective. The crucial necessary condition for a geometric 
optimization problem to be NP-hard is that the solution depends non-locally on 
the data. 

With the LMT-skeleton heuristic, it became feasible to compute minimum-weight 
triangulations fast enough that one could carry out experiments and play with 
various point sets. An instance of a non-local effect was hence discovered by Jack 
Snoeyink |Beirouti and Snoeyink 1998 : the so-called wire is a symmetric polygon 
that does not have a symmetric minimum-weight triangulation, and hence it has 
(at least) two different minimum-weight triangulations, see Figure [2j 






L', ;R 



L', ;R 



Fig. 2. A wire-piece and its two optimal triangulations. The two optimal triangulations are 
labeled as L and R, depending on whether the left (L) or right (R) arms of the isosceles terminal 
triangles, as shown dotted in the right part, belong to the triangulation. This convention is used 
throughout the paper. 
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Fig. 3. A longer wire with several optimal triangulations. Besides the two triangulations of type 
L and R, there are mixed forms where the transition between L and R can occur in intermediate 
positions, as marked by an arrow in the lower example. If the wire is regarded as a point set, 
the shaded pockets outside the wire also have to be triangulated. For each pocket, there are four 
optimal triangulations. However, the pockets are independent of each other, and do not affect the 
properties of the wire. 



Wires can be extended to any length, see Figure [3j A small change in the point 
set at one end of the wire will cause the optimum triangulation to topple globally 
throughout the whole wire. Wires can thus be used to "transmit information" from 
one area of the plane to a remote area (hence the name "wire"). We use vari- 
ous forms of wires as building blocks for constructing gadgets in our NP-hardness 
reduction. 

Such a non-local effect is usually a strong indication that the problem is NP-hard. 
Still, it took almost a decade until another crucial building block was designed: a 
gadget that allows several wires to meet and that carries out "logic" with the 
information that is transmitted by them (the so-called C-connection, see Figure PFT| . 
The design of this gadget heavily depended on computer assistance. We used an 
implementation of the LMT-skeleton heuristic using Otfried Cheong's program Ipe 
ISchwarzkopf 1995| graphical user-interface. 

1.3.2 Computer Assistance. It seems inevitable that a NP-hardness proof for 
the MWT problem requires some amount of computer calculation, since computing 
the weight of a triangulation involves Euclidean distances, and thus square roots. 
Therefore, when comparing two triangulations that differ in more than a pair of 
edges, it is hard to compare their weight by just looking at them. 

However, since the conference version of this paper [Mulzer and R ote 2006J , we 
have made an effort to reduce the part of the proof that depends on computer 
assistance for verification. We have also made other simplifications of the proof 
and made it more accessible, (a) We have tried to simplify and specialize the 
problem as much as possible already at the logical (discrete) level. In particular, the 
PLANAR 1-IN-3-SAT problem is more specialized than in [Mulzer and Rote~2 006j 
and does not need negations, (b) We designed the gadgets in such a way that 
the /3-skeleton edges already form large connected components. This part of the 
construction is trivial to check by computer, and is open even to visual inspection. 
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(c) In the proof, the gadgets split into modules ("pieces") that can be analyzed 
separately, (d) The analysis of each module boils down to checking a small number 
of possibilities, each amounting to computing the minimum- weight triangulation of 
a simple polygon. This can be programmed from scratch in less than an hour's work 
by an experienced programmer, (e) To make the computer-assisted part even safer, 
we have rerun the final calculations with interval arithmetic, using the open-source 
high-level programming language Python, which automatically supports integer 
arithmetic of unbounded length. 

The final computer runs to verify our gadgets took about 10 hours on a relatively 
slow computer. The code is given in Appendix [B] and amounts to approximately 
1000 lines. 



2. OVERVIEW OF THE REDUCTION 

We begin with a high-level overview of the geometric construction. We build our 
point set from a small number of tube-like pieces that fit together at the openings, 
where they share terminal triangles, see Figure HI The boundary edges of the pieces 
must belong to every optimal triangulation, regardless of how we put the pieces 
together (Proposition 16.41 below). If we draw these edges we get a system of tubes, 
separated by polygonal holes. The holes can be triangulated in polynomial time by 
dynamic programming; The tubes will form the essential part of the gadgets that 
model the logical structure of a POSITIVE PLANAR 1-IN-3-SAT formula, whose 
satisfiability is NP-complete to decide, see Section El 

For the analysis, we will show that at each terminal triangle, we have a choice 
of two edges to insert (Lemma 15. ip . When one of the edges is inserted for every 
terminal triangle, the edges enclose a simple polygon, whose optimum solution can 
be computed by dynamic programming. Every piece can therefore be analyzed in 
isolation, by considering a small number of possibilities. 

3. POSITIVE PLANAR 1-IN-3-SAT 

In this section, we describe the POSITIVE PLANAR 1-IN-3-SAT problem, which 
we will use for our reduction. 

Definition 3.1. Let $ be a Boolean formula in 3-CNF. The associated graph of 
<£, G(<£), has one vertex v x for each variable a; in $ and one vertex vq for each 
clause C in <3?. There is an edge between a variable- vertex v x and a clause- vertex 
Vc if and only if x or -ix appears in C. 

The Boolean formula <5 is called planar if its associated graph is planar. 



Lichtenstein [1982] showed that 3-SAT remains NP-complete if the input is re- 



stricted to a planar formula (the PLANAR 3-SAT problem). As Knuth and Raghunathan 1992| 
observed, Lichtenstein's proof implies that it suffices to consider formulae whose as- 
sociated graph can be embedded such that the variables are arranged on a straight 
line, with three-legged clauses above and below them. The edges between the vari- 
ables and the clauses are embedded in a rectilinear fashion (see Figure [5]). 

In our reduction we will use a variant of PLANAR 3-SAT in which there are only 
positive variables and we ask for an assignment to the variables such that in each 
clause exactly one variable is set to true. 
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Fig. 4. A schematic view of the construction: From the five types of building blocks shown at the 
top, one can build a network of tubes like in the lower part. The area within this network is the 
interesting part where the simulation of the logic takes place. The holes inside the structure and 
the pockets between the structure and the convex hull (drawn shaded) form simple polygons that 
can be optimally triangulated in polynomial time. 




Fig. 5. A rectilinear embedding of graph that is associated with the Boolean formula 

(xi V -1x3 V xs) A (-ixi V12V xz) A (x2 V14V -*x*,). 
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Definition 3.2. In the POSITIVE PLANAR 1-IN-3-SAT problem, we are given 
a collection <E> of clauses containing exactly three variables together with a planar 
embedding of the associated graph G(<&) as described above. 

The problem is to decide whether there exists an assignment of truth values to 
the variables of $ such that exactly one variable in each clause is true. 

The POSITIVE 1-IN-3-SAT problem without the planarity restriction was shown 
to be NP-complete by Schaefer 19781. This problem can also be interpreted as a 



SET PARTITIONING problem, where every element of the ground set has precisely 
three sets by which it can be covered. For completeness, we include a proof that 
the planar version of the problem is NP-complete. 

Proposition 3.3. POSITIVE PLANAR 1-IN-3-SAT is NP-complete. 

Proof. It is easy to verify in polynomial time that a given embedding of a 
formula is rectilinear (the validity of an instance) and that a given assignment has 
the l-IN-3 property (the validity of a certificate). Hence, the problem is in NP. 
To show completeness, we describe a reduction from PLANAR 3-SAT. Let I be 
an instance of PLANAR 3-SAT, i.e. a 3-CNF formula <3? and a planar, rectilinear 
embedding of the associated graph. We describe how to transform / into an instance 
of POSITIVE PLANAR 1-IN-3-SAT while maintaining the rectilinear embedding. 
We consider the clauses of one by one. 

If clause C contains only one literal, it can easily be eliminated. If clause C 
contains two literals, say x and y, we can replace it by two three-variable clauses 
(iVi/Vz) and (x V y V -*z), where z is a new variable. 

Before we consider clauses with three variables, we first discuss two useful gadgets 
which enforce equality and inequality between variables in terms of l-IN-3 clauses. 



(a) 



x 



(b) 



x y x 



Fig. 6. The gadgets to enforce inequality (a) and equality (b) between two variables. 

The inequality gadget. The inequality gadget serves to enforce the constraint 
x 44> —iy for two variables x and y. We abbreviate it as (x ^ y). It is implemented 
as 

x y^y = (x,a,y) A (a, b, c) A (a, c, d) A (b, c, d), (1) 

where a, b, c, d are new variables used only inside the gadget (see Figure [6^). Here, 
as always, we denote a l-IN-3 clause by a triple {v\,V2,vz). The clause (^1,^2,^3) 
is satisfied if exactly one of the variables v\, t>2, V3 is true. 
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Lemma 3.4. Given x and y, the expression §\§ is satisfiable iff exactly one of x 
and y is true. 

Proof. The last three clauses enforce that a = 0, because otherwise we get 
b = c = d = 0, and the last clause is not satisfied. The first clause now ensures that 
exactly one of x and y is true. □ 

The equality gadget. Using two copies of the inequality gadget and an extra 
variable a, we can build a gadget to enforce the constraint x 4=> y (see Figure [6b). 
We abbreviate it as (x = y). 



Xi x 2 x 3 x 4 



Fig. 7. Eliminating negations: The variable x in this example has eight connections, some of 
which are negated, as indicated by crosses. We replace i by a chain of variables x\,X2, ■ ■■ with 
alternating truth values and place the connections to the clauses containing x correspondingly. 



Elimination of negated variables. To eliminate the negated variables in the orig- 
inal 3-SAT clauses, we replace a variable x by a chain of variables xi,X2, • • • and 
use the negation gadget to enforce the appropriate relations between them (see 
Figure Ej) . 

Transformation of disjunctions. It remains to handle a disjunctive clause C = 
(x V y V z) with three literals x, y, and z. C can be replaced by 

(x, u, a) A (y, u, b) A (a, b, q) A (u = c) A (d / z) A (c, d, r) . (2) 

See Figure [HJ Note that there is space to accommodate the middle legs of the 
equality and inequality gadgets. 

Lemma 3.5. Given x, y, and z, the expression ([2]) is satisfiable iffx\/y\/z holds. 

Proof. The clause (c, d, r) is equivalent to -icV ->d, since the variable r appears 
nowhere else. Hence c — > ~^d <^=^ u — * z by the third and fourth constraint. Thus, 
satisfiability of the last three constraints is equivalent to u — > z. 

Now, if u = 1, the first three clauses can be satisfied only by x = y = 0, and we 
must have z = 1. The third clause is satisfied by setting a = b = and q = 1. For 
u = 0, the first two clauses reduce to x ^ a and y ^ b. Since q appears nowhere 
else, the clause (a, 6, q) is equivalent to ->a V -*b <^=^ x V y. The value of z can be 
arbitrary in this case. □ 

To make sure that x and z remain reachable from other clauses, we also add two 
additional variables x' and z' and equality constraints. Clauses that were above 
the variables and nested between x and y can now connect to x' instead of x, and 
those nested between y and z can connect to z' instead of z. 
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I I I 
x ■ ■ ■ y ■ ■ ■ z 



\rrTT 



x u q a 

U I 



x 



ir~n 
z c r d 



Fig. 8. The modification for a clause iV y V z with three literals. 

This transformation can be carried out in polynomial time, and by construction, 
the transformed positive formula has a rectilinear embedding, and it has a l-IN-3 
assignment if and only if the original formula is satisfiable. 

4. GEOMETRIC VALIDATION TOOLS 

We summarize the geometric and computational tools we need for establishing 
properties of the gadgets and their minimum-weight triangulations. 

4.1 The /3-Skeleton 

The j3 -skeleton of a point set S is defined as the set of all edges pq between two 
points of S such that the two circles of diameter (3 ■ \pq\ passing through p and q 
are empty, see Figure (9) It provides a sufficient condition for including edges in the 
minimum-weight triangulation: 

Theorem 4.1. The (3-skeleton is a subgraph of the MWT for 



1 + VV27< 1-17682. 



□ 



Theorem 14.11 was first proven bylKeil [19941 for f3 < V2 and was later improved 



by Cheng and Xu [1996 1. The value for (3 is nearly optimal, since there is a lower 



bound of ^5/4+^/1/108 ps 1.16027 |Wang and Yang 2001 1 . In Proposition E2 
we will use the /3-skeleton to identify the boundary edges of our gadgets, which 
must belong to every optimal triangulation. Since the /3-skeleton is defined by a 
local condition, it is very easy to verify that a certain edge is a boundary edge. 
Our gadgets have many points, hence we used a computer to check the /3-skeleton 
property of our boundary edges. The source-code of the straightforward program 
can be found in Appendix IB.3I 

4.2 Triangulating a Polygon by Dynamic Programming 

As we explained above, the analysis of our reduction essentially boils down to 
computing the minimum-weight triangulation of a small number of simple polygons. 



Gilbert [1979 and Klincsek [1980 independently showed how to do this: 



Proposition 4.2. The minimum- weight triangulation of a simple polygon can 
be computed in 0(n 3 ) time. □ 

Proposition 14.21 is proved by a straight-forward application of dynamic program- 
ming, and we implemented a variant of the algorithm in Python. (The source 
code of our implementation is given in Appendix IB.2I ) The dynamic programming 
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Fig. 9. An edge pq of the /3-skeleton. The edge ab does not belong the MWT since the diamond 
test is violated. 

algorithm assumes an idealized model in which arithmetic with edge lengths can 
be done exactly. This means that sums of radicals (i. e. Euclidean distances) need 
to be compared in constant time. However, it is not known how to do this even 
in polynomial time [Blomer 1991J. Therefore, our program uses interval arithmetic 
to compare sums of distances, and for all our gadgets it computes an unambiguous 
answer. 

Our program identifies the edges in a minimum-weight triangulation T inside a 
given simple polygon P. We call the total weight of these edges the internal cost 
of P and denote it by c(T). The internal costs computed by our program will be 
discussed extensively in Section 16. 1L It is easy to modify the program so that it 
computes the optimum under the restriction that certain edges are forbidden. This 
is needed for the proof of Lemma 15.11 

5. TERMINAL TRIANGLES 

The point set is constructed from small elementary pieces that fit together at ter- 
minal triangles. These terminal triangles are isosceles triangles xyz with base yz, 
that lie symmetrically with respect to a coordinate axis, and they come in two 
sizes (small and large). The coordinates of vertical terminal triangles are given in 
Table HI See Figure [10] or [16HTT1 for an illustration. These triangles can be rotated 





X 


y 


z 


5 


small 
large 


(0.00, 0.00) 
(0.00, 0.00) 


( -2.7 , 11.2 ) 
(-11.61, 48.16) 


( 2.7 , 11.2 ) 
(11.61, 48.16) 


<5 sm all = 5.655172 

^large = 24.06 



Table I. Terminal triangles. The "difference terms" S will described below in Section 15. II 



by multiples of 90°, and translated by multiples of 0.01 in each coordinate. The 
large triangle is a scaled copy of the small triangle, by a factor of 4.3. 

The basic pieces that we construct are point sets with two or three terminal 
triangles. The /3-skeleton edges will form a simple polygon through these points, 
with a missing edge at each terminal triangle. We will show that in each terminal 
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Wieft 



(2.7,11.2)/ 



r 




Wright 



V, 



a 



Fig. 10. The point set W in Lemma 15.11 



triangle, the edge xy or the edge xz must be present in every minimum-weight tri- 
angulation (Lemma 15. ip . Thus, in order to determine how a piece with k terminals 
can be triangulated in a minimum-weight triangulation, it suffices to consider 2 fc 
possibilities for the positions of each terminal edge. For each possibility, we find 
the optimum solution by dynamic programming. 

Let W be the point set given in Figure [TUl which forms part of the wire mentioned 
in the introduction (see Section [1.3.ip , It is symmetric about the y-axis. The points 
vo, x, and v' lie on the x-axis, and the points u and v! lie 0.1 units below this line. 
The left and right halves of the lower boundary are equal, and they are equal to 
the upper boundary (turned by 180°). We denote by Wi e f t and Height the left half 
and the right half of the set W, up to and including the points x, y, and z. 

Lemma 5.1. Let P be a point set that contains a (translated, scaled, rotated) 
copy of the point set W, but no other points in the shaded polygon of FiaureMCK Any 
minimum-weight triangulation T of P that contains all edges shown in Figure Q3 
must contain at least one of the edges xy or xz. 

Proof. The proof distinguishes 21 cases and deals with each case by a small 
calculation. 

Since in any triangulation, the point u must be incident to some edge that em- 
anates into the lower half-plane, T must contain at least one of the edges uv\, UV2, 
. . . , uvq. Similarly, T must contain at least one of the edges u'v[, u'v' 2 , ■ ■ ■ , u'v' 6 . 
Now, for each i,j = l,...,6, the edges uvi and u'v 1 - enclose a simple polygon. For 
each polygon, we calculated, by dynamic programming, the optimum triangulation 
as well as the optimum triangulation that uses none of the edges xy and xz. By 
symmetry, we have to consider only 1 < i < j < 6. The results are shown in 
Table IE □ 

We remark that the lemma is somewhat robust against perturbations and remains 
true even if all points of W are moved by a small distance. Figure [11] shows 
that small perturbations do not make triangulations feasible or infeasible. The 
triangulations in Table [XT] contain at most 18 internal edges. If each point is moved 
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Vl,Vj 

4.00304 




V2,v 2 
4.67187 



V3,V 3 
4.00308 



4.67191 



3.83677 



4.67195 



Vi,V 2 
4.50551 



Vl,V 3 
4.00306 



Vi,V 4 
4.50553 



Vl,V B 
4.00308 

Vl,v' 6 
4.50555 

4.50553 

4.67189 

v 2 ,v' s 
4.50555 

4.67191 



[202.72577, 
202.72594] 




181.98535, 
181.98550] 




[159.65909, 
159.65922] 




137.62231, 
137.62242] 



[116.25971, 
,116.25980] 



[94.68883, 
> 94.68890] 




192.35556, 
192.35572] 




[181.19243, 
181.19258] 




170.17404, 
170.17418] 




159.49274, 
159.49287] 

148.70730, 
148.70742] 

170.82222, 
170.82236] 

159.80383, 
159.80396] 

149.12253, 
149.12265] 

138.33709, 
138.33720] 




[198.72255, 
198.72273] 




[177.31332, 
177.31348] 




[155.65587, 
155.65601] 




[132.95028, 
132.95040] 



[112.42284, 
> 112.42294] 



[90.01680, 
90.01688] 




[187.84988 
187.85005' 




[177.18921, 
177.18937] 




[165.66836, 
165.66851] 




155.48952, 
155.48966] 

144.20162, 
144.20175] 

[166.31654, 
166.31669] 

[155.13180, 
155.13194] 

144.61685, 
144.61698] 

[133.66506, 
133.66518] 




[198.72255, 
198.72273] 




155.65587, 
155.65601] 



[112.42284, 
> 112.42294] 




177.18921, 
177.18937] 



"3>«4 
4.50555 


xfS- [148.64070, 
Z\Z/\ 148.64082] 


Xf7\ [144.13502, 
L^L^jK 144.13515] 


4.00310 


X~^\ [137.95940, 
/\/\ 137.95951] 


/r=X [133.956X8, 
/V)/\ 133.95630] 


4.50557 


xf^V [127.17396, 
/\rT> 127.17406] 


/=T\ [122.66828, 
122.66839] 


4.50557 


,<f^V [126.94101, 
f^VV-A 126.94111] 


/T=X [122.43533, 
/VV\ 122.43544] 


V4,l/ e 
4.67193 


[116.15557, 
L^)rT> 116.15566] 


/V7\ [111.48354, 
C^~XT> 111-48364] 


«s,«4 

4.50559 


,<^V [105.47427, 
(\rT> 105.47435] 


^=TX [100.96859, 
100-96868] 



Table II. The optimum triangulations with and without the edges xy and xz, for each polygon 
bounded by edges uvi and u'v'^, as indicated in the left column. The second column shows an 
optimal triangulation without the edges xy and xz, together with the interval for its internal cost, 
calculated to 5 digits after the decimal point. (There are always at least two optimal triangulations, 
since the upper diagonal in the middle can be nipped. The program has picked one solution whose 
lower interval bound is smallest.) The remaining columns list all candidates for optimal solutions. 
One can check that all these solutions contain the edges xy and xz. (In fact, one can see that, 
if there are two optimal solutions in a row, they must have exactly the same weight since they 
contain parts which can be replaced by symmetric parts.) The lower bound on the difference 
between the optimum in the third column and the optimum in the second column is given in the 
first column. The smallest difference (3.83677) occurs for the case v$,v' 5 . The numbers in this 
table were computed with the PYTHON program in Appendix IB. 41 
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Fig. 11. The point set W in Lemma 15. II is in sufficiently general position. 

by at most a := 0.04, the length of each edge changes by at most 2a, and the weight 
of any triangulation changes by at most 36a = 1.44. 

Since the difference between the optimal triangulation and the best triangulation 
without the edges xy and xz is bigger than 3.83677 > 2 • 1.44, the claim of the 
lemma is not invalidated if each point is moved by at most 0.04 in any direction. 
We will use this observation in Section [61 where we describe a version of the wire- 
piece in which the offset between the terminal triangles is slightly larger than in 
the standard wire-piece. This makes it possible to build wires of arbitrary length 
(see Lemma l6.ip . 

5.1 Putting pieces together 

When putting pieces together, we make sure that Lemma 15.11 can be applied at 
the junction: at every terminal triangle, the piece contains a copy of half of the 
point set W (possibly rotated, and possibly scaled by 4.3). This allows us to split 
the input set into components at the terminal triangles. For each terminal triangle 
xyz, we need to consider just the two choices xy and yz for the terminal edges. We 
denote these choices by L and R, depending on whether we use the left or right 
arms, as viewed from the tip of the terminal triangle. 

When determining the behavior of a piece with k terminals in a minimum-weight 
triangulation, we have to compare the optimal triangulations for the 2 k choices for 
the positions of the k terminal edges. However, these triangulations cover different 
areas and cannot be compared directly. Depending on whether the terminal edge 
xy or the terminal edge xz is chosen at some given boundary, the triangulated 
area excludes or includes the area of the terminal triangle itself, and the optimal 
triangulation can thus be expected to be cheaper or more expensive. To offset this 
effect, we define the reduced weight of a triangulation T, c(T), by adding a penalty 
or subtracting a bonus term from the actual weight. This will allow us to see 
directly which configurations of terminal edges are better than others, and which 
configurations cannot possibly be part of an optimal triangulation. 

More precisely, we define a difference term S that depends on the size of the 
triangle (large or small) , see Table HI When comparing the costs of different tri- 
angulations of a piece or of different pieces, we only look at the internal cost, the 



Technical Report B05-23 (revised), arXiv:cs/0601002v3, February 2008 



Minimum-Weight Triangulation is NP-hard ■ 15 




Fig. 12. All elementary pieces 

cost of the internal edges, since the boundary edges are fixed. When the triangu- 
lated area of a piece includes the terminal triangle, we subtract 5 from the actual 
weight; otherwise we add S. The precise values of <5 sma n and <5i arge are not important 
and have been chosen to make the proofs convenient. The following lemma holds 
irrespective of the values of <5 S maii and <5i arge - 

Lemma 5.2. Consider a layout of pieces S where every terminal triangle is 
shared by two pieces. For each terminal triangle, fix one of its two long sides. 
Consider the minimum-weight triangulation T of S that is restricted to contain 
those fixed edges as well as all piece boundaries. For each piece, calculate the re- 
duced internal weight of the minimum-weight triangulation of the piece bounded by 
the given fixed edges. 

Then the overall reduced internal weight of all pieces differs from the total weight 
of the triangulation T by a constant that is independent of the choice of the fixed 
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thickening adapter 




Fig. 13. Enlarged view of the small elementary pieces 

edges. 

Proof. In total, the effect of adding and subtracting 5 cancels on the two sides 
of each terminal triangle. The cost of the terminal edge (xy or xz) itself is not 
accounted for in the sum of internal costs, but since xy and xz have the same 
length, this amounts to a constant difference. The boundaries are also fixed, and 
so is the cost of triangulating the pockets and the holes. □ 

We will see below (Proposition 16. 4j ) that all boundary edges are part of the MWT, 
since they belong to the /3-skeleton. We know from Lemma 15. II that, by choosing a 
terminal edge from each terminal triangle in all possible ways, we are guaranteed 
to find the optimal triangulation. The consequence of Lemma 15.21 is that we need 
only look at the reduced (internal) cost when comparing these choices. 
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Fig. 14. Enlarged view of the central part of the C and Co connections. 

6. ELEMENTARY PIECES 

We have 10 elementary pieces, shown in Figures [T2]414t The (i) wire-piece (Fig- 
ure [16]) is the main tool achieving a non-locality for conveying information over 
long distance. As mentioned in the introduction, such wire-pieces were originally 
conceived by Jack Snoeyink [ Beirouti and Snoeyink 1998| . In order to construct 
wires of arbitrary length, we also have (ii) an extended wire-piece. (Figure [TBI) . The 
horizontal offset between the two terminal triangles of a wire-piece is 27.4, i.e., 
with these wire-pieces alone, one can bridge distances which are multiples of 27.4. 
The offset for the extended wire-piece is 82.21 = 3 • 27.4 + 0.01. Thus, with the 
right combination of wire-pieces and extended wire-pieces, one can form a straight 
connection of arbitrary length, provided it is long enough: 

Lemma 6.1. Consider two small vertical terminal triangles at the same height 
with a horizontal distance d > 230000 that is a multiple o/O.Ol. Then the two 
triangles can be connected by a sequence of wire-pieces and extended wire-pieces. 

An analogous statement holds for vertical connections. 

Proof. If the distance between two terminal triangles is d = 0.01 • z for some 
integer z > 3 • 2740 • 2739 = 22 514 580, they can be connected by concatenating 
y := z mod 2740 extended wire-pieces with [^/2740j — 3y wire-pieces. □ 

The work-horse of our gadgets is (hi) the C-connection and its mirror image, 
(iv) the C -connection. These are the only pieces with three terminal triangles, two 
large ones and a small one. They serve two purposes: they allow us to introduce 
branches into the network of wires, and they effect "negation" of the information 
that is transmitted through the wires. We also have a (v) Co-connection with only 
two terminals. It serves as a placeholder for the C or C connection when the third 
terminal is not needed. 

The remaining pieces are variations of the wire-piece for building more flexible 
wire shapes: The (vi) left bend and the (vii) right bend allow us to introduce 
arbitrary turns into wires. The vertical part of the left bend corresponds to a wire- 
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Fig. 15. Tracing a wire through a rectangular corridor, forming a so-called C-C'-link (see Sec- 
tion E!) The dotted parts are long enough to allow connections of any length (Lemma I6.1|l . 

piece that is rotated 90° counter-clockwise, whereas the vertical part of a right bend 
corresponds to a wire-piece that is rotated 90° clockwise. We need both a left and 
a right bend, because each piece fits to a given terminal triangle of another piece 
only at one end. 

With these pieces, we can arrange wires to follow any rectangular layout, provided 
we blow the layout up sufficiently, see Figure HSJ 

Lemma 6.2. Given a C-connection piece and a C -connection piece and a rect- 
angular path between them, the small triangles of the two pieces can be connected 
by a sequence of wire-pieces, extended wire pieces, left and right bends, within the 
corridor of width 2000 around the path, provided that the path contains a straight 
portion of length at least 250 000, both in the horizontal and in the vertical direction, 
and the corridor does not interfere with other pieces or intersect itself. 

Proof. The left and right bend both fit into a box of size 700 x 700, and by 
Lemma 16-H we can bridge any distance larger than 230 000 in the horizontal and 
vertical direction. The wires can be positioned sufficiently far away from the bound- 
ary of the rectangular strip to make sure that the boundaries of different connections 
do not interfere with each other. □ 

We also have (viii) a thick left bend, which is just a scaled copy of the left bend. 

Finally, we have two resizing wire-pieces: the (ix) thickening adapter and its 
mirror image, the (x) thinning adapter. They are used to interpolate between small 
and large terminal triangles in the clause gadgets. 

Figures [16] and [T7] show the coordinates of the wire-piece and the extended wire- 
piece. The complete data for the pieces are available on the InternetE] 

Proposition 6.3. The elementary pieces have the following properties: 
All coordinates are multiples of 0.0001. Each piece contains two or three copies 
of a terminal triangle, together with a copy of W\ e f t or Wright {Figure[W§ , possibly 
rotated by 90°, and possibly scaled by 4.3. The coordinates of the terminal triangles 
are multiples o/O.Ol. □ 

Each piece has two terminal triangles, with the exception of the connection pieces, 
which have three terminal triangles, as shown in the figures. The edges shown in 



1 http : //www. inf . fu-berl in. de/inst/ag-ti/people/rote/Sof tware/MWT /all-coordinates. txt 
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Fig. 16. The coordinates of a wire-piece. 




Fig. 17. The coordinates of the extended wire-piece. 

the figure, with the exception of the two equal sides of each terminal triangle, form 
the boundary of the pieces. 

Proposition 6.4. The boundary edges of all elementary pieces are contained in 
the (3-skeleton, for (3 = 1.1806. This remains true when the pieces are connected 
using their terminal triangles, as long as the boundary part of different pieces are 
sufficiently far apart. Hence, all boundary edges belong to the MWT. 

Proof. The lemma can be verified manually by inspecting the relevant regions 
for the boundary edges. Since the number of boundary edges is very large, we have 
implemented a computer program which computes the /3-skeleton. The source code 
can be found in Appendix IB.3I □ 
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The boundary decomposes into two or three connected boundary pieces. As can 
be seen in Figure [Til some boundary pieces are not just paths, but they form trees. 

6.1 Analysis of the Pieces 

In this section we discuss the exact properties of our gadgets as they were computed 
by our computer program. We present tables which show the internal cost, the 
reduced internal cost c (see Section 15. ip and the relative reduced internal cost c 
(i.e., the difference between the reduced cost of a certain configuration and the 
minimum reduced cost incurred by any configuration of a gadget). 




Fig. 18. Optimal solutions for all cases for the wire-piece and the extended wire-piece 



For example, let us look at the extended wire-piece. It has two terminal triangles. 
For each of the four combinations of states, LL, LR, RL, and RR, the right part 
of Figure [THl shows an optimal solution, together with the internal cost. Table ITTT1 
summarizes this information. The second column lists the number of optimal so- 
lutions (or potential optimal solutions, i.e., all solutions which the program could 
not exclude from being an optimal solution, with the given precision of the calcu- 
lation). Indeed, looking back at Figure [THl one can see that the solution for case 
LR has a symmetric copy with the same cost. The astute reader may wonder why 
the multiplicities for cases LR and RL are not larger. For example, it looks as if 
the isosceles u xyz" triangle for case LR could be placed also at the very left end 
or at the very right end. However, this is not the case, since the three repetitions 
of the wire-piece that constitute the extended wire-piece are not exact repetitions: 
the ends have been moved apart, stretching the middle piece and thus causing a 
deviation from the regular periodical pattern. The third column of Table IIIII lists 
the internal cost, repeating the information from Figure [THl 



pattern 


multiplicity 


internal cost c 


reduced internal cost c 


relative reduced cost c 


LL 


1 


455.471523 435 


455.471523435 


0.000 000000 


LR 


2 


466.990 265 006 


455.679 921006 


0.208 397 570 


RL 


1 


444.283 180 745 


455.593 524 745 


0.122 001310 


RR 


1 


455.471523 435 


455.471523435 


0.000 000000 



Table III. Analysis of the extended wire-piece 
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The fourth column is the reduced internal cost c, obtained by adding or sub- 
tracting 5 sm aii, as appropriate. For example, for case LR, we have to subtract 
25 S maii = 11.310 344, since this triangulation covers the larger area on both the left 
and the right end. For cases LL and RR, addition and subtraction of 5 sma ii cancel, 
and the reduced cost remains unchanged. 

The rightmost column reduces all numbers by subtracting the column minimum, 
to make the differences more visible. It is now plain to see that LL and RR have 
the same cost, which is also obvious by symmetry. We denote these relative reduced 
costs by c, and we think of them as penalties for deviating from the "ground state" 
where the smallest reduced cost is achieved. They are non-negative by definition. 

The corresponding table for all pieces, together with pictures of the optimal trian- 
gulations, is given in Table IVTIl of Appendix [A] Table [Tvl condenses the information 
for all pieces with two terminals. 

Since the pieces are symmetric, LL and RR always have the same cost. The only 
exceptions are the resizing pieces. Let E\ ~ 0.000051 be the relative reduced cost 
of the thickening adapter in state LL. 





wire-piece 


extended wire-piece 


thickening adapter 


thinning adapter 


LL 


0.000 000 000 


0.000 000 000 


0.000 051402. . . = ei 


0.000 000 000 


LR 


0.210 506 663 


0.208 397 570 


0.018 887 246 


0.018 887 246 


RL 


0.125 593 246 


0.122 001310 


0.014627 250 


0.014 627 250 


RR 


0.000 000 000 


0.000 000 000 


0.000 000 000 


0.000 051402. . . = Ei 





Co connection 


left bend 


right bend 


thick left bend 


LL 
LR 
RL 
RR 


0.000 000 000 
0.044001701 
0.020 571757 
0.000 000 000 


0.000 000 000 
0.086 460 895 
0.125 593 246 
0.000 000 000 


0.000 000 000 
0.210 506 663 
0.125 593 246 
0.000 000 000 


0.000 000 000 
0.891261046 
0.020 571757 
0.000 000 000 



Table IV. The relative reduced costs c for the pieces with two terminals 



The C-connection (and its mirror image, C) is the only piece with three termi- 
nals. The results for this piece is shown in Table [V] and in more visual form in 
Figure [19] We encode the states by three letters ABc, where the capital letters A 
and B refer to the left and right large terminal triangles (L or R), as usual, and c 
refers to the small terminal triangle (1 or r). 



c 


internal cost c 


reduced internal cost c 


relative reduced cost c 


a 


LL1 


14 027.752 986 494 


14 033.408158 494 


0.003 861076. . 


■ = 5i 


RRr 


LLr 


14039.059 470 993 


14 033.404298 993 


0.000 001575. . 


■ = E2 


RR1 


LR1 


14075.894120134 


14 033.434292134 


0.029 994 716 




LRr 


LRr 


14087.200604 634 


14 033.430432 634 


0.026135 215 




LR1 


RL1 


13 979.654697176 


14 033.424869176 


0.020 571 757 




RLr 


RLr 


13 990.965 042 750 


14 033.424870 750 


0.020 573 332 




RL1 


RR1 


14027.749125 419 


14 033.404297419 


0.000 000 000 = 





LLr 


RRr 


14039.064100 296 


14 033.408 928 296 


0.004630 878. . 


. = s 2 


LL1 



Table V. Analysis of the pieces C and C . Since C is the mirror image of C, it is represented in 
the same table 
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c 


relative reduced cost 


a 


LL1 




^? 


0.003 861076.. 


■ = 6i 








RRr 


LLr 






0.000 001575.. 


■ = £2 




RR1 


LR1 




0.029 994 716 


A. 






LRr 


LRr 






0.026135 215 




LR1 


RL1 


R zC 


^? 


0.020 571757 


R z 






RLr 


RLr 




x^ 


0.020 573 332 


R z 


fey 


> 


RL1 


RR1 




0.000 000 000 = 





\ 


fey 


y 


LLr 


RRr 


R zC 


x^ 


0.004 630 878.. 


. = s 2 




fey 


y 


LL1 



Fig. 19. All cases for the analysis of the pieces C and C . 

All tables were computed using fixed-precision interval arithmetic with 15 decimal 
digits after the decimal point (i.e., multiples of 10 -15 ). The resulting intervals 
were rounded to 9 digits for displaying. It turned out that the intervals were 
small enough so that they could be rounded to unique 9-digit numbers. Thus, for 
example, 0.000 051402 denotes a number that is guaranteed to lie strictly between 
0.000 051 401 5 and 0.000 051 402 5. 

By inspecting Tables ITVl and PVT one can see that there is a penalty for having two 
triangles of different states. If both terminal triangles of a piece with two terminal 
triangles, or the two large terminal triangles of C or C, are in a different state, we 
call this a breach. A breach is so expensive that it will never occur in an optimal 
solution. We summarize this information for later use. 

Proposition 6.5. In every piece, a breach has relative reduced cost at least 63 := 
0.01. If there is no breach, the relative reduced cost is zero, except for C, and C , 
and the resizing wire-pieces. For the resizing wire-pieces, the relative reduced cost 
is at most e x < 0.000051. □ 

It follows that long wires have a preference to be triangulated uniformly: 

Lemma 6.6 (The Wire Lemma). Let T\, T2 be two small terminal triangles 
connected by an arbitrary sequence of wire-pieces, extended wire-pieces, and left 
and right bends. If T\ and T2 are in the same state, all the connecting pieces are 
in that state and the relative reduced cost is 0. If T\ and T2 are in different states, 
the relative reduced cost is at least 83 = 0.01. □ 

The states of C and C without a breach are LL1, LLr, RR1, and RRr. From 
Table El one sees that the preferred states are LLr, and RR1, where the small 
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triangle has the opposite state from the large triangles. We call these two states 
the consistent states, and the states LL1 and RRr inconsistent. We define 

e 2 := c(C, LLr) = c(C, LLr) - c(C, RR1) = c(C, LLr) - c(C, RR1) - 2<5 smal i 

8 1 := c(C, LL1) = c(C, LL1) - c(C, RR1) = c(C, LL1) - c(C, RR1) 

8 2 := c(C, RRr) = c(C, RRr) - c(C, RR1) = c(C, RRr) - c(C, RR1) - 2<5 smal i, 

where c(C, LLr), c(C, LLr), etc. denote the cost of the C-connection in the respective 
configuration, see Table IVl 

Proposition 6.7. In C or C , every consistent state has relative reduced cost 
at most E2 < 0.000 002. The two inconsistent states have cost 8i ~ 0.003 861 and 
82 ~ 0.004 631. Thus, if there is no breach in C or C , the relative reduced cost is 
at most 82- □ 

In general, we denote by £1, 62, "small" quantities that we would like to neglect, 
and which can be made arbitrarily small by refining the coordinates of the pieces 
or the digits of <5 sma n and <5i arge . On the other hand, Si, 82, ■ ■ ■ , denote "large" 
quantities whose difference will be made productive for the proof. 

7. LARGER GADGETS 

We now describe how to assemble the elementary pieces from the last section to 
obtain larger gadgets which model the logical structure of a PLANAR 1-IN-3-SAT 
formula. We need variables, clauses, and the connections between them. The main 
building block for the variables and clauses is the bit loop (Section IT. X [) . Each bit 
loop stores a logical state (L or R). To ensure consistency between different bit 
loops, we use C-C -links (Section I7.2p . By connecting several bit loops, we build 
clauses (Section IT.3f) and variables (Section 17. 4D . The optimal triangulation of a 
clause is obtained if and only if exactly one of its inputs is in state L (the l-IN-3 
property). A variable is modeled as a chain of bit loops connected by C-C"-links 
so that all loops are in the same state. Variables can be connected to clauses using 
the wire-pieces from the last section. 

7.1 The Bit Loop 

The bit loop is used to represent a logical state (L or R). It if formed by connecting 
any selection of four connection pieces (C, C or Co) into a loop, using four thick 
bends, as shown schematically in Figure [201 

Lemma 7.1. In an optimal triangulation, the 8 large terminal triangles of a bit 
loop are in the same state (L or R). 

Proof. The only interaction with other pieces is through the at most four small 
"exit" triangles. Let us keep the states of these triangles fixed. If we set all 8 
large terminal triangles to L or to R, the maximum relative reduced cost is at most 
4 • 8 2 < 0.02 = 25 3 , see Proposition ED 

On the other hand, if the 8 large terminal triangles are not all equal, there are at 
least two breaches. Hence, by Proposition 16.51 the relative reduced cost is at least 
283, and such a triangulation cannot be optimal. □ 
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C" 




C 



Fig. 20. An instance of a bit loop with two C-connections and one C'-connection. At the left C- 
connection, a right bend is attached, as the beginning of a wire joining it to another C'-connection. 
The drawing is to scale, but the shapes are simplified, 

We can therefore simply refer to the state of the bit loop as the common state of 
its large terminal triangles. Generally, when it is clear that there is no breach, the 
state of a C, C or Co-connection refers to the state of its large terminal triangles. 

7.2 The C-C'-Link 

A C-connection can be joined to a C'-connection by an arbitrary sequence of wire- 
pieces, extended wire-pieces, and left and right bends (Lemma 16. 2p . see Figure [T5l 
We call the two connections at the end together with their joining wire a C -C -link. 

C-C-links are used in two contexts: The most common case is when the C- 
connection of the C-C'-link belongs to a bit loop V\ and the C'-connection belongs 
to another bit loop V-2- In the clause gadget (Section IT.3[) . we also have C-C- 
links where the C-connection is connected to a thickening adapter and a thinning 
adapter, while the C connection is part of a bit loop. By Lemma 17.11 we know 
that in each of V\ and V?, the large triangles are in a uniform state (L or R), and 
in Section T7.3I we will see that an analogous statement holds for the large terminal 
triangles of the thickening and the thinning adapter. We thus have to consider four 
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cases for these states, as shown in Table IVIi For each combination, we can work 
out the optimum state of the small terminal triangles by considering Table |V1 (The 
first two and last two rows of Table [V] are sufficient.) The optimal choice for the 
thin wire is shown in the third column. 



Vi 


v 2 


wire 


state of C 


state of C 


optimum relative reduced cost 


L 


L 


r 


LLr 


LLr 


c(C,LLr) + £(C",LLr) = e 2 


L 


R 


1 (or r) 


LL1 


RR1 


c(C, LLr) + c(C", RRr) = e 2 + <5i 


R 


L 


1 (or r) 


RR1 


LL1 


£(C,RRr) +c(C", LL1) = 8 2 


R 


R 


1 


RR1 


RR1 


c(C, RR1) + £(C", RR1) = e 2 



Table VI. The best state of the wire for a C-C'-link. For each given combination of states (first 
two columns), the last column gives the resulting cost. 

Lemma 7.2. Each C-C'-link incurs a relative reduced cost of at least £2- The 
cost reaches this minimum if and only if the large terminal triangles in the two pieces 
have the same state. Otherwise, the cost increases to 5\ + £2 or 62, respectively, see 
TableWl 

Proof. We need only exclude the possibility that the two small terminal trian- 
gles are in different states. By the Wire Lemma (Lemma I6.6p . this would cause a 
relative reduced cost of at least £3 = 0.01, which is larger than any of the costs in 
Table ED □ 

We call the C-C'-link inconsistent if the large triangles in the two connection pieces 
have different states. An inconsistent C-C'-link incurs a relative reduced cost bigger 
than 8\. 

7.3 The Clause Gadget 

The clause gadget is formed by 3 pairs of bit loops, a and a, (3 and (3, 7 and 7, 
see Figure [2TJ The two loops of each pair are connected by a C-C'-link in which a 
thickening adapter, a C-connection piece labeled DOWN, and a thinning adapter 
is interspersed, and another similar C-C'-link with a Co-connection piece in the 
middle. The a-a group has three external connections: DOWN, UP, and ENTRY. 

The DOWN C-connection piece of the pair a-a is connected to the UP C con- 
nection of (3, and so on in a circular way. Finally, each of a, (3, and 7, has a 
C-connection piece, labeled ENTRY, which will be connected to one of the "input" 
variables of the clause by a C-C'-link. 

Lemma 7.3. In an optimal triangulation, a and a are in the same state {L or 
R), and the two large terminal triangles in the DOWN C-connection are in the 
opposite state. 

Analogous statements hold for (3,(3 and 7,7. 

Proof. By symmetry, it suffices to consider a and a. The a-a group has three 
external triangles UP, DOWN, and ENTRY which connect it to the outside world. 

There are two "ground states" with no internal breaches or inconsistencies. In 
these states, a and a (including UP and ENTRY) have the same state and DOWN 
and the upper Co-connection piece in the middle are in the opposite state. 
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Fig. 21. Schematic view of the clause gadget. The dashed parts of the wires are long enough for 
sufficiently many copies of the wire-piece and the extended wire-piece to ensure that the wire can 
reach from the lowest C connection to the left C connection of & (Lemma 16. l|l . 

If we ignore the three external connections, the relative reduced cost of these 
two states is 2ex + 2e2- The contribution from the two C-C'-links is 2^2- The 
contribution from the four resizing wire-pieces is 2e\ since we have two thickening 
and two thinning adapters: the "penalty state" for the thickening adapter (LL) is 
the optimal state for the thinning adapter, and vice versa. 

Now, fix a setting for the three small external triangles. If we select the ground 
state which is consistent with the majority of the three exit triangles, we have 
at most one inconsistent exit connection, causing an additional penalty of at most 
52+2^2 (two consistent connections with at most £2, and one inconsistent connection 
with at most 62)- 

Thus, there is always a ground state solution with cost at most 82 + 2e\ + 4^2 < 
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0.005. If a and a are not in the same state, or if DOWN is in the same state 
as a and a, there must be at least two inconsistencies (or even breaches) in the 
connections between a, a, and DOWN, causing a cost of at least 28\ > 0.007, which 
is bigger than for the ground state solution. □ 

Since DOWN is always in a different state from UP, ideally, to successive loops 
in the cyclic sequence a(3j should be in a different state. However, since there 
are three loops, there has to be at least one inconsistency, which provides us the 
asymmetry necessary for a l-IN-3-clause. 

Thus, we obtain: 

Lemma 7.4. The reduced cost of the clause gadget (excluding the three ENTRY 
C -connections) achieves its minimum if and only if exactly one of a, (3, and 7 is in 
state L. Any other triangulation incurs a cost that is at least 84 := 0.0007 larger. 

Proof. By Lemma [7,31 it suffices to analyze the 8 possible configurations of the 
pairs a-a, /3-/3, and 7-7 (ignoring the contribution of the ENTRY connections to 
the clause). 

The "internal" contribution or each pair of 2e2 + 2ei from its resizing pieces and 
its two internal C-C"dinks is constant, and thus we can ignore this amount when 
comparing the various possibilities. 

Let us look at the C-C'dink between two successive pairs, say a-a and /3-/3. If 
they are equal, they cause an inconsistency, and the relative reduced cost is 82 if 
a = (3 = L and 5\ + £2 if a = /3 = R, according to Table IVII If a and (3 are in 
different states, the link is consistent, and the cost is £2- 

Since the situation is unchanged under cyclic shifts of the sequence a/37, it is 
enough to consider four cases: 

— a/37 = LLL or a/37 = RRR: In these cases, we have three inconsistencies, and 

the relative reduced cost is bigger than 3<5i > 0.01. 
— a/37 = LLR: the relative reduced cost is 82 + 2e 2 > 0.0046. 
— a/37 = RRL: the relative reduced cost is 81 + 3e2 < 0.0039. 

□ 

In the previous lemma, we have ignored the relative reduced cost of the ENTRY 
C-connections. They will be accounted for as part of the C-C'dinks that they form 
with the variables, to be described next. 



4 C' C k C 4 C 




Co ▼ C C T C C t 



Fig. 22. Schematic view of a variable gadget with three upper and three lower exits. The pattern 
can be repeated arbitrarily. 
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7.4 The Variable Gadget 

The gadget that represents a variable of the l-IN-3 instance is formed by chaining 
together sufficiently many bit loops and connecting them in a row using C-C-links, 
as shown in Figure 1221 The dotted arrows denote potential exit wires that establish 
C-C'-links with the corresponding clause gadgets of the previous section. Unused 
exits can be replaced by Co connections. 

Lemma 7.5. Consider a variable chain V together with the C-C'-links that con- 
nect it to the clauses, including the ENTRY C -connections in the clauses. The 
reduced cost of this point set achieves its minimum if and only if all the bit loops 
and the C -connections in the clauses are in the same state. Any other triangulation 
incurs a cost that is at least 5i larger. 

Proof. This follows directly from Lemmas 17.11 and 17.21 □ 

8. THE REDUCTION 

Theorem 8.1. Minimum-weight triangulation is strongly NP-hard. 

Proof. A rectilinear embedding of the given PLANAR 1-IN-3-SAT formula can 
be constructed on a grid of size 0{n) x O(n). The reduction procedure then simply 
replaces the edges, variables and clauses of the PLANAR 1-IN-3-SAT formula by 
the appropriate gadgets, taking care to leave enough space between the individual 
wires. 

This procedure yields a point set S. By construction, the boundaries of all 
the gadgets are part of the /3-skeleton of S (Proposition 16. 4p . and therefore they 
belong to the minimum-weight triangulation. The faces outside the wires are simple 
polygons and can be optimally triangulated using dynamic programming. For each 
gadget, we know the desired "ideal triangulation" (Lemmas l6.6l and l7.1H7. 5] ) and can 
calculate its weight. Adding up these weights and the weights of the faces outside 
the wires yields a target weight w. By construction, the input instance is l-IN-3 
satisfiable if and only if the minimum weight of a triangulation of S is w. Otherwise, 
the weight of the shortest triangulation is at least w + 0.0007 (Lemmas 17.41 and IT.5|) . 

The set S is a subset of an 0(n) x 0(n) grid, but it does not fill the whole area: 
it follows the linear structure of the edges of the rectilinear embedding. Therefore, 
S has 0(n) points. Hence the triangulation has 0(n) edges of length 0{n). By 
calculating all edge lengths with an absolute error of 0(l/n 2 ), the reduction algo- 
rithm can thus calculate, in polynomial time, a threshold w such that the input 
formula is satisfiable iff there is a triangulation of length at most w. □ 

9. CONCLUSION 

9.1 Running Times for Computer Verification of the Proof 

In designing our gadgets and our proof, we have tried to achieve a balance between 
the number of different pieces, which affects the complexity of the human-readable 
part of the proof and the number of case distinctions, and the size of the pieces, 
which affects the complexity (running time) for the mechanical part of the proof 
that has to be checked by computer (or accepted by faith). 

The running time is dominated by the 0(n 3 ) dynamic programming algorithm 
for triangulating simple polygons. The largest point sets that we handle have 493 
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points (the left bends). To total time to run all verifications (with exact integer 
interval arithmetic) was about 10 hours on a four-year-old moderate PC. 

9.2 Dimensions. 

Our gadgets have constant size, but they consist of several thousands of points and 
are quite enormous. For example, the clause gadget (see Figure l2Tj) has dimensions 
on the order of 250 000 x 250 000, with coordinates that are specified as multiples of 
10 -4 . On the other hand, the difference between a satisfiable and an unsatisfiable 
SAT instance is reflected in a minute difference of 0.0007 in the MWT cost. 

With some work, it would be possible to reduce this to some more "reasonable" 
figures or to amplify the difference between satisfiable and unsatisfiable instances, 
but we did not find it worth the effort to do so. First of all, the current gadgets 
are already a result of tedious experiments, pushing points into various directions 
and trying to understand what happens. Some parts of the design, in particular 
the C-connections, are very delicate. Secondly, the dimensions would still be very 
large. One can certainly reduce the constant 230 000 of Lemma 16. II by providing a 
greater variety of extended wire-pieces, but the bit loop, for example, (Figure [20]) 
already has size approximately 7000 x 7000, and it does not seem easy to push the 
size very much below these limits, unless one comes up with a completely different 
design. 

9.3 Open Problems 

Several interesting problems remain open. First of all, it is not known whether the 
MWT problem is in NP, since it is not known how to compare sums of Euclidean 
lengths in polynomial time [Blomer JMQ, but this difficulty is more of an algebraic 
nature. To define a variant of MWT which is in NP, one can take the weight of 
an edge e as the rounded value HNhl- With appropriate scaling, our proof also 
establishes NP-completeness for this variant. 

Our reduction shows that it is NP-hard to approximate the MWT with a relative 
approximation error which is better than 0(l/n 2 ): The difference between a satisfi- 
able instance of PLANAR 1-IN-3-SAT and an unsatisfiable instance is reflected in a 
constant increase of the MWT cost, and, as mentioned in the proof of Theorem 18. 11 
the total cost of the MWT is 0(n 2 ). 

One can probably reduce this bound to 0(n log n), and thus establish that it is 
NP-hard to achieve a relative approximation error better than 0(l/(nlogn)), by 
using the fact that the interior of a convex k-gon of perimeter p can be triangulated 
with weight 0(p\ogk). First, the wires and all gadgets form linear structures of 
total length 0(n); thus, the length of the gadget boundaries, and the MWT inside 
the gadgets is only O(n). This leaves the holes to be triangulated. It should be 
quite straightforward to extend our construction in such a way that, apart from 
0{n) constant-size holes, only 0{n) convex holes with a total of 0(n) vertices 
remain: one would insert paths of additional points whose /3-skeleton separates the 
big holes from the jagged wire boundaries and cuts the holes into convex pieces, 
apart from "linear" structures that cover only 0(n) area. 

These non-approximability results do not rule out the existence of a polynomial- 
time approximation scheme. For a long time, attempts to extend techniques from 
geometric approximation algorithms to the MWT problem have only led to constant 
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factor approximations (see Bern and Eppstein 1996 for a survey). Remy and Steger [2006 



showed that it is possible to compute a (1 + e)-approximation of the MWT in time 
n O(iog n) ^ providing strong evidence that a PTAS might exist. 

In practice, the LMT-skeleton heuristic is extremely fast in computing LMTs. 
Combined with with bucketing techniques and fast preprocessing techniques |Drysdale et al. 1995 



one empirically achieves almost linear running times. Thus, in this respect the 
MWT problem seems to be similar to the Knapsack Problem, which is also NP- 
hard but easy to solve in practice [Keller er et al. 2007] . It would be interesting 
to analyze the LMT-skeleton heuristic for random point sets. The good practical 
performance indicates that the expected running time for random inputs might be 
polynomial, or even close to linear. On the other hand, as point sets get huge, they 
will contain, with non-negligible probability, some larger and larger point configu- 
rations that are hard to triangulate. 
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A. THE DATA FOR ALL PIECES 

Table IVlII shows, for all pieces, the internal cost, the reduced cost, and the relative 
reduced cost, for all configurations of boundary triangles. Figures I23rf28l as well as 
Figure [TBI show for each configuration an optimal triangulation. The two resizing 
wire-pieces are symmetric, and thus only one of them, the thickening adapter, is 
given. When applying the results to the thinning adapter, the labels must be 
nipped: LL becomes RR and RR becomes LL. The raw costs for the thick left 
bend are the same as for the left bend, scaled by 4.3. But since the large terminal 
triangles use a different value of 5, the reduced costs bear no direct relation to those 
of the left bend. The last column is what is summarized in Table IIVI 

We extended our Python programs to output their results in the form of tables 
(in TpjX format) and figures. As the output format for the figures, we used the xml- 
format of the Ipe prograrrfl [ Schwarzkopf 1995] , which can be directly converted to 
eps or pdf format or further edited. Many illustrations of this paper were obtained 
in this way. In particular, the table and all figures of this section, as well as 
Figures [TBI and Table ITU were generated directly from our programs without manual 
intervention. 



2 http : //tclab .kaist . ac .kr/ipe/ 
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pattern 


multiplicity 


internal cost c 


reduced internal cost c 


relative reduced cost c 


wire-piece 


LL 


1 


144.135078832 


144.135078832 


0.000000000 


LR 


2 


155.655929495 


144.345585495 


0.210506663 


RL 


1 


132.950328078 


144.260672078 


0.125593246 


RR 


1 


144.135078832 


144.135078832 


0.000000000 


extended wire-piece 


LL 


1 


455.471523435 


455.471523435 


0.000000000 


LR 


2 


466.990265006 


455.679921006 


0.208397570 


RL 


1 


444.283180745 


455.593524745 


0.122001310 


RR 


1 


455.471523435 


455.471523435 


0.000000000 


thickening adapter 


LL 


1 


7188.980062460 


7207.382390460 


0.000051402 


LR 


1 


7237.113898304 


7207.401226304 


0.018887246 


RL 


1 


7177.684294307 


7207.396966307 


0.014627250 


RR 


1 


7225.784667058 


7207.382339058 


0.000000000 


C connection 


LL1 


1 


14027.752986494 


14033.408158494 


0.003861076 


LLr 


1 


14039.059470993 


14033.404298993 


0.000001575 


LR1 


1 


14075.894120134 


14033.434292134 


0.029994716 


LRr 


1 


14087.200604634 


14033.430432634 


0.026135215 


RL1 


1 


13979.654697176 


14033.424869176 


0.020571757 


RLr 


1 


13990.965042750 


14033.424870750 


0.020573332 


RR1 


1 


14027.749125419 


14033.404297419 


0.000000000 


RRr 


1 


14039.064100296 


14033.408928296 


0.004630878 


Co connection 


LL 


1 


12733.864882577 


12733.864882577 


0.000000000 


LR 


2 


12782.023884279 


12733.908884279 


0.044001701 


RL 


2 


12685.770454334 


12733.885454334 


0.020571757 


RR 


1 


12733.864882577 


12733.864882577 


0.000000000 


left bend 


LL 


1 


5425.386907832 


5425.386907832 


0.000000000 


LR 


2 


5436.783712726 


5425.473368726 


0.086460895 


RL 


2 


5414.202157077 


5425.512501077 


0.125593246 


RR 


1 


5425.386907832 


5425.386907832 


0.000000000 


right bend 


LL 


1 


5271.525633165 


5271.525633165 


0.000000000 


LR 


2 


5283.046483828 


5271.736139828 


0.210506663 


RL 


2 


5260.340882410 


5271.651226410 


0.125593246 


RR 


1 


5271.525633165 


5271.525633165 


0.000000000 


thick left bend 


LL 


1 


23329.163703675 


23329.163703675 


0.000000000 


LR 


2 


23378.169964722 


23330.054964722 


0.891261046 


RL 


2 


23281.069275432 


23329.184275432 


0.020571757 


RR 


1 


23329.163703675 


23329.163703675 


0.000000000 



Table VII. Analysis of the pieces 



B. COMPUTER PROGRAMS 

We have written programs that check all of the claimed properties. (Proposi- 
tions l6.3H6.4p . If any condition fails, the programs raise an exception. To give 
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Fig. 23. Optimal solutions for all cases for the Co connection 




Fig. 24. Optimal solutions for all cases for the C connection 

an idea of what was automatically checked, we show an excerpt of the log-file, 
concerning the extended wire-piece (Figures 1171 and [TH]) . 

======================== extended wire-piece ==================== 

All coordinates are multiples of 0.0001 
terminal triangle basepoint: (-41.100,0) 

edge vector in state L: (-2.700,11.2) 
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Fig. 



26. Optimal solutions for all cases for the left bend 
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RR: 5271.525633165 



Fig. 27. Optimal solutions for all cases for the right bend 

edge vector in state R: (2.700,11.2) 
terminal triangle basepoint: (41.110,0) 

edge vector in state L: (-2.700,11.2) 

edge vector in state R: (2.700,11.2) 
All terminal coordinates are multiples of 0.01 
duplicate point (s). 

The point set is symmetric with respect to the vertical axis x=0.005. 

cos(alpha)~2 = 3120343/10000000 = 0.312034; beta = 1.205637. 

case LL: 44 points. 455.471523435 455.471523435 

case LR: 45 points. 466.990265006 455.679921006 

case RL: 43 points. 444.283180745 455.593524745 

case RR: 44 points. 455.471523435 455.471523435 

The complete Python source code, as well as the data for the pieces, are avail- 
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RR: 23329.163703675 



Fig. 28. Optimal solutions for all cases for the thick left bend 

able on the Internet^] In total, there are about 1000 lines of code. The programs 
are mostly rather straightforward, so we only show four modules: The module 
mwt_polygon.py contains the dynamic programming algorithm for computing the 
MWT of a simple polygon (Section IB.2|) . It relies on the module arithmetic .py 
for the representation of fixed-precision decimal quantities as (long) integers, for 
(rudimentary) interval arithmetic with integers, and for calculating the Euclidean 
length as an interval (Section IB. 1|) . The module check_beta.py determines the 
smallest value j3 for which the given boundary edges belong to the /3-skeleton (Sec- 
tion lB.3j ). The module check_W.py checks all cases for the point set W for proving 



3 http : //www. inf . f u-berlin. de/inst/ag-ti/people/rote/Sof tware/MWT/python/ 
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Lemma 15.11 (Section IB.4|) . 
B.l Arithmetic 

Since the dynamic programming algorithm involves only additions and compar- 
isons, interval arithmetic is not a big deal. We use fixed-point arithmetic, which is 
simulated by integer arithmetic. Hence all calculations are exact, and we need not 
care about directed rounding when doing arithmetic with intervals. 

The only non-trivial part is the calculation of Euclidean distances by the length 
function, which involves a square root. We use the decimal package of Python, 
which provides an arbitrary-precision square root operation, to obtain a starting 
approximation of the square root. The approximation is then checked and refined 
directly, in a straightforward way. 

from decimal import * 
cdigits = 5 

scale = 100000 # initial values 

def setprecision(computedigits) : 
global cdigits , scale 
cdigits = computedigits 
scale = 10**computedigits 

def length (pl,p2): 

"""Euclidean distance between two points 
pi and p2 are points with integer coordinates. 
The results is an integer interval.""" 
l_squared = (p2 [0] -pi [0] ) **2 + (p2 [1] -pi [1] ) **2 

# obtain an initial approximation using the 

# square root operation of the decimal package 
oldprec = getcontextO .prec 

getcontextO .prec = len(str (l_squared) )/2+2 

a = int (Decimal (l_squared) . sqrtO .quantize(Decimal(l) , 

rounding=RDUND_FL0DR) ) 
getcontextO .prec = oldprec # restore the old precision 

# correct and verify the result . 

# The final result will be correct regardless of how good 

# the initial approximation was (although it may take longer) . 
while a**2 > l_squared: a -= 1 

while (a+l)**2 <= l_squared: a += 1 

if a*2 < l_squared: return Interval_integer (a, a+1) 

else: return Interval_integer(a,a) 

def decimal_to_scaled_int (s) : 
return int(s*scale) 

def print_scaled_int (n,ndigits=None) : 

"""print integer with decimal point inserted 
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"ndigits" digits from the right end""" 
if n<0: 

sign = "-" 

n = -n 
else : 

sign = ' ' 

if ndigits==None : ndigits=cdigits # default 

return sign + C7„d.°/„0"+str(ndigits) + "d") °/„ divmod(n, 10**ndigits) 

class Interval_integer(object) : 

Him represents an interval [a,b] with a<=b """ 

def init (self , a,b=None) : 

self . a = a 
if b==None : 

self.b = a 
else : 

self.b = b 

def add (self,y): 

if isinstance (y , (int.long)): 

return Interval_integer (self .a + y, self.b + y) 
else: # assume y is an Interval_integer : 

return Interval_integer (self .a + y.a, self.b + y.b) 

def quantize(self , ndigits) : 

"""try to find a unique value which is guaranteed 

to be the rounded value of the interval""" 

scalehalf = 10**ndigits/2 

rounded_value = (self . a/scalehalf + l)/2 

val = rounded_value*2*scalehalf 

if val-scalehalf < self .a and val+scalehalf > self.b: 

return rounded_value 
print self, "cannot be rounded uniquely to", 
print ndigits, "fewer digits." 
raise ValueError 

def quantize_and_print (self .ndigits , initialdigits=None) : 
"""print the unique value which is the rounded value 
of the interval with "ndigits" remaining digits, if it exists""" 
if initialdigits==None : initialdigits=cdigits # default 
return print_scaled_int ( 

self . quantize (initialdigits-ndigits) , ndigits) 

def str (self) : 

return " ["+str(self . a)+" , "+str (self .b)+"] " 
def repr (self) : 

return " Interval ("+repr (self .a)+" , "+repr(self .b)+") " 

def keep_best (1) : 

"""keeps a list that is guaranteed to contain the smallest 
value from a list of intervals""" 
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threshold = min(x.b for x in 1) 

return [ x for x in 1 if x.a <= threshold ] 



B.2 Dynamic Programming for Triangulating a Polygon 

The procedure mwt of the module mwt_polygon.py implements the classical 0(n 3 ) 
dynamic programming algorithm for optimally triangulating a simple polygon p 
[Gilbert 1979} iKlincsek 19 80], It has an optional argument excluded for specifying 
edges that cannot be used in the triangulation. The input polygon p is assumed to 
have integer coordinates. Euclidean edge lengths are calculated as integer intervals 
using arithmetic . length () . 



In contrast to |Klincsek [19 80] , our procedure does not explicitly test the edges 



for crossings; it only tests whether all triangles that are used in the triangulation 
are oriented counterclockwise. If the input is a simple polygon, this is sufficient 
to ensure that the resulting triangulation is non-crossing, see Lemma IB. II below. 
(Hence, the program will also triangulate certain non-simple polygons, but it will 
not triangulate polygons which are oriented clockwise!) 

To show that the our constructions do not depend on assumptions about han- 
dling point sets which are not in general position, triangles with three points on 
a line are (temporarily) considered as valid triangles of a triangulation, but these 
triangles, as well as any (partial) triangulations that contain such triangles, are 
flagged as degenerate. If such a triangulation would "survive" as a candidate for 
an optimal solution, the checking routines would report an error (see for example 
the procedure check_W_cases () in Section rB,4|) . Triangles with coinciding vertices 
are not considered. 

Lemma B.l. Let T be a set of triplets (i,j,k) with l<i<j<k<n such 
that for a convex polygon P = p\ . . .p n , the triangles PiPjPk for (i,j, k) £ T form a 
triangulation of P. 

Let P be an arbitrary simple polygon. If all triangles PiPjPk for (i,j,k) G T are 
oriented counter-clockwise, they form a triangulation of P. 

Proof. This can be seen by counting the number of triangles in which a given 
point x of the plane is contained. This number can only change when x crosses 
an edge of a triangle. However, all triangle edges have another triangle on the 
opposite side, with the exception of the triangle edges that are edges of P. Thus, 
the number of triangles in which a point of the plane is contained is constant except 
at the boundary of P, where it changes by ±1. Since this number is when the 
point x is far away, every point x in the interior of P is covered by exactly one 
triangle, and no triangle sticks out of P. □ 

import arithmetic 

def mwt(p,excluded=[] ) : 

"""compute minimum-weight triangulation of a simple polygon p 
by dynamic programming, using integer interval arithmetic. 
The coordinates of p are integers. 



Some edges may be excluded from consideration. 
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The polygon p must be oriented counter-clockwise! 
Otherwise, no solution will be found. 
The result is a list of triplets 

(weight, solution, degenerate_f lag) 
where weight is given as an interval, sorted by the 
lower bound weight. a of the solution quality 
To show that the result does not depend on assumptions 
about handling point sets which are not in general position, 
degenerate solutions (triangles with three points on a line) 
are considered, but they are flagged as degenerate.""" 

opt = {} ## array of optimal (partial) solutions. 

n = len(p) 

for i in range(n-l): # initialization 
o = arithmetic. Interval_integer (0) 

## calculate the weight only of interior edges 
o.info = None # backtracking information and flag is stored 

# within the Interval_integer objects 
o .degenerate_f lag = False 
opt[i,i+l] = [o] 

# opt[i,j] stores a LIST of all potential candidate solutions. 

for 1 in range(2,n): # l="length" of the edges 
for i in range (n-1) : 
j = i+1 

# check for triangulations containing edge (i,j): 
if (i,j) in excluded or (j,i) in excluded: 

## we are interested in MWT that does not contain those edges 

continue 

if p[i]==p[j]: # this point appears twice on the boundary 

continue 
if (i,j) == (0,n-l): 

edgelength = # final edge is not counted 
else : 

edgelength = arithmetic . length(p [i] ,p [j] ) 
for k in range (i+1, j) : 

if opt .has_key ( (i ,k) ) and opt .has_key( (k, j) ) : 
ori = orientation (p [i] ,p [k] ,p [j] ) 
if ori >= 0: 
inl = opt[i,k] 
in2 = opt [k , j ] 

allcand = [ computesol(inl,in2,nl,n2,k,ori==0, edgelength) 

for nl in range (len(inl)) 

for n2 in range (len(in2)) ] 
opt[i,j] = arithmetic. keep_best (opt . get((i, j) , [] )+allcand) 

# sort final result by lower bounds of intervals: 
allopt = [(x.a,x) for x in opt [0, n-1]] 

allopt . sort () 

opt [0, n-1] = [x[l] for x in allopt] 

# retrace optimal solution: 

return [ (arithmetic . Interval_integer(x. a, x.b) , 

retrace (opt, 0, n-1, ind) , x.degenerate_f lag) 
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for x,ind in zip( opt [0,n-l] , range (len(opt [0,n-l] ) ) ) ] 

def orientation(pl ,p2 ,p3) : 
x2 = p2[0]-pl[0] 
y2 = p2[l]-pl[l] 
x3 = p3[0]-pl[0] 
y3 = p3[l]-pl[l] 
return x2*y3-x3*y2 

def computesol(ll,12,il,i2,k,flag,edgelength) : 
r = 11 [il] + 12 [i2] + edgelength 
r.info = (k,il,i2) # k=splitting point. 

# il,i2: pointers into sublists 
r .degenerate_f lag = (flag or 11 [il] .degenerate_f lag 

or 12 [i2] .degenerate_f lag) 

return r 

def retrace (opt, i, j ,ind) : 

"""recursively recover the solution, based on the "info" 

fields stored with the solution values in the array "opt". 
Recover the solution for position "ind" in the list "opt". 
The result is a list of triplets (i,k,j), each representing 
the two edges p[i]p[k] and p[k]p[j].""" 
if opt [i ,j] [ind] . info : 

k,il,i2 = opt [i , j] [ind] . info 

return [ ( i , k , j ) ] +retrace (opt , i , k , i 1 ) +retrace (opt , k , j , i2) 
else : 

return [] 

B.3 Checking the /3-skeleton 

This procedure is straightforward: for every edge pq of the boundary, it runs 
through all remaining points r and checks whether they violate the /3-skeleton 
condition, by calculating the (squared) cosine of the angle a = prq with the cosine 
law. The running time is 0(n) per edge of the boundary, thus at most 0(n 2 ) in 
total. The relation between the angle a and the ratio (3 between the diameter of 
the circumcircle of pqr and the distance \pq\ is given by 

sina = l//3, or cos 2 a = 1 — 1//3 2 

To understand this program, one has to know how the pieces are given. Each 
piece is represented as a dictionary. The field piece ["partlist"] is a sequence 
of (names of) boundary parts, usually just ("lower", "upper"); only the C- 
connection has three parts. Then piece ["lower"] and piece ["upper"] (or what- 
ever the names in the "partlist" are) are the actual boundary parts, as lists of 
points. Each point is a pair of precise Decimal numbers as provided by the decimal 
package. 

from math import sqrt 

from checking_routines import smallest_unit 
from basic_routines import map_coordinates 
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beta_threshold =1.18 # ">=1. 17683 !" 
betaprecision = 10**7 

def coordinate_to_int (c) : 

"scale coordinate to integer" 

return int (c/smallest_unit) # smallest_unit=0 . 0001 

def length2(p) : 

"""squared length of vector ' 
return p[0]**2 + p[l]**2 

def checkbeta(piece) : 

min_cos_squared = 100*betaprecision # infinity 
for part in piece ["partlist"] : 

# scaling all coordinates does not affect beta-skeleton: 
poly = map_coordinates (coordinate_to_int , piece [part]) 
for i in range (len (poly) -1) : # run through all edges: 
pi = poly[i] 
p2 = poly[i+l] 

# check edge (pl,p2): 

# run through all points p: 

for part_2 in piece ["partlist"] : 

for p in map_coordinates (coordinate_to_int .piece [part_2] ) : 
if p<>pl and p<>p2 : 

# compute cos_squared, representing the 

# squared cosine of the angle alpha=pl ,p,p2 
xl = (pl[0]-p[0], pl[l]-p[l]) 

x2 = (p2[0]-p[0], p2[l]-p[l]) 
innerproduct = xl[0]*x2[0] + xl[l]*x2[l] 

# so far, all calculations were exact, 
if innerproduct>0 : 

cos_squared = ( innerproduct**2 * betaprecision / 

(Iength2(xl)*length2(x2)) ) 

# scaled by betaprecision and rounded down; 

# yields a conservative (lower) estimate of beta 
min_cos_squared = min(min_cos_squared, cos_squared) 

min_cos_squared_unscaled = min_cos_squared/f loat (betaprecision) 
beta = 1/sqrt (l-min_cos_squared_unscaled) 

# This final calculation is just in double-precision floating-point, 
print "cos(alpha)~2 = '/.d/'/.d = %8.6f; beta = */.8.6f." 7, ( 

min_cos_squared, betaprecision, min_cos_squared_unscaled, beta) 
if beta < beta_threshold: 

print "beta too small: ", beta, ">=1. 17683 !", beta_threshold 

raise RuntimeError 
return beta 

B.4 Proving Lemma [5711 

This program runs through all 21 subpolygons of the polygon W in the proof of 
Lemma 15.11 For illustration, we give here the output produced for the first case: 

Case vl, vl': difference = 

Technical Report B05-23 (revised), arXiv:cs/0601002v3, February 2008 



44 



W. Mulzer and G. Rote 



4 . 00304 

Best solution value without terminal edges: 

[202.72577,202.72594] 
There is/are 2 best solution(s) without restriction: 

[198.72255,198.72273] 

[198.72255,198.72273] 

There is an extended version (not shown here) which also generates a drawing 
of the triangulations, and which has been used to generate Table HU The module 
solve_all_patterns .py, by which the data of Appendix lAl were generated, is 
similar. 

# checks presence of one ot two "terminal" edges in thin wire 

from basic_routines import map_coordinates 
import arithmetic 

from arithmetic import print_scaled_int , decimal_to_scaled_int 
import mwt_polygon 

from pieces import W_upper, W_lower, W_point_x, W_point_y, W_point_z 

# W_upper and W_lower form the set W. 

def check_W_cases () : 

mindiff = 999 * arithmetic . scale # infinity 

for i in ranged, 7): 

for j in ranged, 7): # for K=i<=j<=6: 

p = W_lower[i:len(W_lower)-j] + W_upper 
i_x = p . index (W_point_x) 
i_y = p.index(W_point_y) 
i_z = p . index (W_point_z) 

excl = [d_x, i_y) , (i_x, i_z)] # critical edges 

pcalc = map_coordinates (decimal_to_scaled_int , p) 

# First, compute optimum solution without restriction 
elist = mwt_polygon.mwt (pcalc) 

for (result , sol , degenerate) in elist: 

# convert list of triplets to list 

# of pairs representing edges: 

sole = [(x,y) for [x,y,z] in sol] + \ 

[(y,z) for [x,y,z] in sol] 
if not ( (i_x, i_y) in sole or (i_x,i_z) in sole or 
(i_y,i_x) in sole or (i_z,i_x) in sole): 
print "incorrect : " ,i, j , result, sol, degenerate 
print ("optimal solution does not contain" , 

i_x,i_y, "or", i_x,i_z) 
raise StandardError 
if degenerate: 

print "DEGENERATE SOLUTION :", i , j , result , sol .degenerate 
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raise StandardError 

# Now, compute optimum solution with the two edges excluded: 
elist2 = mwt_polygon.mwt(pcalc, excluded=excl) 

# consider just the first solution from the list: 

diff = elist2 [0] [0] . a - max(result.b for (result, s,f) in elist) 
mindiff = min(mindif f ,dif f ) 

print "Case v'/.d, v'/.d' : difference =" 7. (i,j) 
print print_scaled_int (dif f ) 

print " "*10+"Best solution value without terminal edges:" 

print_scaled_interval (elist2 [0] [0] ) 

print " There is/are" , len(elist) , 

print "best solution(s) without restriction:" 

for result, sol, f in elist: print_scaled_interval (result) 

print 

print "Minimum overall difference =", print_scaled_int (mindiff ) 

def print_scaled_interval(x) : 
print " "*14+" [7.s,7„s] "7„(print_scaled_int (x . a) , print_scaled_int(x.b)) 

print "Program check_W.py. Check terminal edges xy and xz in W" 
arithmetic . setprecision(5) # 5 digits are sufficient. 
check_W_cases() 
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